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Cn ' Abstract. Many applied time-dependent problems are characterized by an 

f^.,^' additive representation of the problem operator. Additive schemes are con- 

rrt ■ structed using such a splitting and associated with the transition to a new 

time level on the basis of the solution of more simple problems for the indi- 
vidual operators in the additive decomposition. We consider a new class of 
additive schemes for problems with additive representation of the operator at 
the time derivative. In this paper we construct and study the vector operator- 
difference schemes, which are characterized by a transition from one initial the 
evolution equation to a system of such equations. 



(N 



> 

00 



< 

Q I Introduction 

For the approximate solution of multidimensional unsteady problems of mathe- 
matical physics there are widely used different classes of additive schemes (splitting 
schemes) fSlitSlTT] . Beginning with the pioneering works [21,16] the most simple way to 
construct additive schemes is in the splitting of the problem operator on the sum 
f~^ , of two operators with a more simple structure — alternating direction methods, 

CSJ ' factorized schemes, predictor-corrector schemes etc. [12] . 

jy-s , In the more general case of multicomponent splitting, classes of uncondition- 

f^ ' ally stable operator-difference schemes are based on the concept of summarized 

^^ , approximation. In this way, we can construct the classic locally one-dimensional 

schemes (componentwise splitting schemes) 018], additively-averaged locally one- 
dimensional schemes [3l[12]. 

A new class of unconditionally stable schemes — vector additive schemes (mul- 
ticomponent alternating direction method schemes) is actively developed (see, eg, 
C^ ' [lEl] ) ■ They belong to a class of full approximation schemes — each intermediate 

problem approximates the original one. The most simple additive full approxi- 
mation schemes are based on the principle of regularization of operator-difference 
schemes. Improving the quality of operator-difference schemes is achieved using 
additive or multiplicative perturbations of operators of the scheme [7]. Regular- 
ized additive schemes for evolutionary equations of the first and second order are 
constructed for equations as well as systems of equations [T31IIS]- Both the stan- 
dard schemes of splitting with respect to separate directions (locally-onedimensional 
schemes), splitting with respect to physical processes and regionally-additive schemes 
based on domain decomposition for constructing parallel algorithms for transient 
problems of mathematical physics jHfTOlIIl] . 
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At present, different classes of additive operator-difference schemes for evolution- 
ary equations are constructed via additive splitting of the main operator (connected 
with the solution) onto several terms. For a number of applications it is interesting 
to consider problems in which the additive representation demonstrates an operator 
at the time derivative. In this work, for this new class of evolutionary problems 
the vector additive operator-difference schemes are constructed and studied. The 
work is organized as follows. Section 1 provides a statement of the problem along 
with a simple a priori estimate of the stability for the solutions with respect to 
initial data and right-hand side. This estimate is nothing but our reference point 
when considering the vector problem and the operator-difference schemes. The 
vector differential problem is considered in Section 2. The central part of the work 
(Section 3) deals with the construction and investigation of the stability of vector 
additive schemes. Possible generalizations of the results are discussed in Section 4. 

1. Statement of the problem 

Let _ff be a finite-dimensional Hilbert space, and A, B, D be linear operators in 
H. We consider grid functions y of finite-dimensional real Hilbert space H, for the 
scalar product and norm in which we use the notations: (•, •), \\y\\ = {y, y)^'^. For 
D — D* > we introduce space Hd with scalar product {y,w)D — {Dy,w) and 
norm ||y||D = {Dy,y)^/'^. 

In the Cauchy problem for evolutionary equation of first order we search function 
y{t) £ H, which satisfies the equation 

fill 

(1.1) B^ + Au = f{t), i>0 

and the initial condition 

(1.2) u(0) = u° 

at given f{t) G H. 

We assume that linear operators A and B, acting from H into H {A : H ^ H, 
B : H ^ H), are positive, self-adjoint and stationary, that is 

For problem p. II) . (jl.21) we can obtain different a priori estimates, which express 
the stability of the solution with respect to the initial data and right hand side in 
different spaces. We restrict ourselves to the simplest of them, trying to get the 
same type of estimates for both the scalar and vector problems as well as for the 
solution of both differential and difference problems. 

Multiplying scalarly both sides of equation (|l.ip in H by u, we get 

2 df *'^"' "^ ''' ^^^' ^^ ^ ^•^' "'*■ 
For the right hand side we use the estimate 

if,u)<iAu,u) + ^{A-'fj). 

This yields the following a priori estimate for the solution of problem (jl.ip , (|1.2p : 

t 

(1-3) \\uml<hYB + lJ\\fis)\\Uds, 
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which expresses the stabihty of the solution with respect to the initial data and 
right hand side. 

Standard additive difference schemes are characterized by decomposition (split- 
ting) of the operator A onto the sum of operators of a simpler structure. For 
example, we assume that for operator A we have the following additive representa- 
tion: 

p 

(1.4) A = ^^,, A„^AI>0, a = 1,2,..., p. 

Q = l 

Additive difference schemes are based on the basis of (|1.4I) . where the problem is 
decomposed into p subproblems. The transition from time level t" to the next level 
t"-+i = i" -I- T, where r > is the time step and y" = y(i"), t" = nr, n = 0, 1, ..., 
is associated with solving problems for individual operators Aa,a = l,2,...,p in 
additive decomposition (|1.4p . 

The subject of our consideration will be another case. In a number of problems 
the computational complexity is not associated with operator A, but with operator 
B at the derivatives in time. In this case, to decrease the computational complexity 
of problem (jl.ip , (|1.2p we employ the additive representation 

p 

(1.5) B = J2Bc., B^ = B*^>0, a = 1,2,..., p. 

a=l 

instead of (|1.4p . The transition to a new time level is connected with the solution 
of some auxiliary Cauchy problems for equations 

Ba,-^+Aua, = Uit), t>0 a = l,2,...,p 
with specified appropriate initial conditions. 

2. Vector problem 

By definition, put u = {ui, U2, ..., Wp}. Each individual component is defined as 
the solution of similar problems 

(2.1) £i?;3^-(-Aw„ = /(t), t>0, 

(2.2) w,(0)=w°, a = 1,2,..., p. 

Here is the simplest coordinate-wise estimate for the stability of the solution. 
Subtracting one equation from another, we get 

A{ua~Ua^i) =Q, a = 2, 3,..., p. 

Taking into account the positivity of operator A this gives 

Ua = Ua-l, a = 2,3, ...,p. 

For separate component u^ we obtain the same equation as for u: 

Y,Bp'^+Auo. = f{t), i>0, a = 1,2,..., p. 
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For the same reason, there are a priori estimates 

t 
(2.3) Wu^ml < hYs + lj\\fm\-^ds, a = l,2,...,p. 



It foUows that 

Ua{t) = u{t), t > 0, a — 1,2,..., p. 

Therefore, as the solution of original problem (ll.ip . (|1.2p we can take any component 
of the vector u(i). 

For the vector evolutionary problem we can obtain a priori estimates for vector 
u, considering the problem in Hilbert space H = H^ with the scalar product 

p 

(U,v) = ^{Ua,Va). 
a=l 

This technique is used, for example, in |12) when considering additive schemes with 
splitting pTi)) . 

We rewrite equations (|2.ip in the form 



dujs 

/3 = 1 



B„^-i^S^^+S„u„ = /„(i), i>0, a = l,2,...,p. 



where /„ = Bo,A ^ f . This allows us to write the system of equations in vector 
form 

(2.4) C^+Du==f. 

Operator matrix C and D have the form 

(2.5) C = {a^}, C^p = B^A~^Bp, 

'D = {Dafi}, Da/j ^ Ba6ap, a, /3 = 1, 2, ..., p, 

where Sa^ is the Kronecker delta. Equation (|2.4p is supplemented by the initial 
condition 

(2.6) u(0) = u°. 

The principal advantage of notation (|2.4p results from the fact that 
C = C* > 0, D = D* > 

inH. 

Here is a priori estimate for the solution of vector problem (|2.4I) - (|2.6I) . This 
estimate, on the one hand, is more complicated than (|2.3p and, on the other hand, 
we will use it as the guideline in the consideration of the operator-difference schemes. 

Multiplying both sides of (|2.4p scalarly in H by du/dt, we get 

f ^du du\ 1 d ,„ , f^-du 

^'■'^ (^^'d^)+2d^(^"'")n''^ 

Taking into account (|2.5p . we obtain 

^ \ /3=1 /3=1 
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and for the right hand side of (|2.7p we have 

Similarly ([O]), ^^, from (HJ]), dlH) it follows the estimate 



(2.9) l|u||^<||u"l|2, + iy'||/(.)l|^_.d.. 



Taking into account (|2.5p . we have 



p 

a=l 



Thus, estimate (|2.9p can be considered along with (|2.3p as the vector analogue of 
estimate (|1.3p . Taking into account p.Sp , estimate (|2.7p gives the stability of any 
individual component of vector u(i). 

3. Additive vector schemes 

Splitting schemes for the approximate solution of (jl.ip . (|1.2p . ()1.5p will be con- 
structed on the basis of usual schemes with weights for vector problem (|2.ip , (|2.2p . 
The standard two- level scheme with weights for problem (jl.ip , (|1.2p has the form 

(3.1) B^ y_ + A{ay^+^ + [l - a)y^) = ip^ , n = 0,l,..., 



r 
where, for example, 

^» = /(at"+i + (1 - a)r), 

and a is a weight parameter (usually < cr < 1). 

In the general theory of operator-difference schemes stability developed by A. A. 
Samarskii [SHlOj. there were obtained the exact (unimproved) stability criteria for 
two-level and three-level operator-difference schemes in various norms. They can 
be directly used in the study of schemes with weights (|3.ip . Here is a typical result. 

Theorem 3.1. If cr > 1/2, then operator- difference scheme 113. 1\) is absolutely 
stable in Hb and for the difference solution the level-wise estimate is valid 

(3.2) ll2/"+'lll<l|y"lll + p¥'"ll^-i. 

Proof. By definition, put 

1 / 1 \ 7;"+l 

ff(n) _ ^,,"+1 I n _ rr\-,,^ — _^,"+i ' ~-"^ . - I - ^ \ y 



y^in, ^ ^yn+, ^ (^ _ ^^yu ^ _^yn+, + ^n) ^ ^ |^^ _ _ 

Multiplying scalarly in H both sides of (|3.ip by y'^(") , we get 

^(i?(y"+i-2/"),y"+i+2/") + 

r (a - i) ^B y''^'-y\ ^!1^^ + (Ay-("),y'^(«)) = (^",y^W). 
For the right hand side we use the estimate 
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If cr > 1/2, we obtain desired estimate p.2|) for the stability of the numerical 
solution with respect to the initial data and right hand side, which is the grid 
analog of estimate (|1.3p for the solution of problem (|l.ll) . (jl.2[) . This concludes the 
proof. D 

To solve vector problem (|2.ip . (|2.2p we apply the following difference scheme: 

/ yn+l _ n n _ n-1 

(3.3) B^ie^ ^ + (1-61)^^^ — ^^^ 



P y-Ti _ n-1 

n = 0,l,..., a == l,2,...,p. 

Unlike p. 11) ) scheme (l3.3p is a three-level one and has two weight factors and a. 

Numerical implementation of scheme (j3.3p is associated with sequential solving 
grid problems 

{eBl + aTA)yl+^^xl. a^l,2,...,p 
with transition from time level t" to new time level f"+^. For vector additive 
scheme p.3p it is possible to implement a parallel organization of computations — 
an independent calculation of the individual components. 

Using notation (|2.5p . we write operator-difference scheme p.3p in the vector form 

n+l _ 2v" + V"~^ 

(3.4) 0G^ ^—^ + 



n _ ^n-1 

C^^ i + D(f7y"+i + (1 - 2a)y" + ay"-i) 



where 

G = {Gap}, Gafi ~ BaA Ba6aj3, 

g" = {52}, 5S = Ba^-V", a,/3 = l,2,...,p. 
Thus, in ([0|) operator G = G* > 0. 
Taking into account that 

yTC _ y"— 1 y"+i _ y"^i y"+i — 2v"^"'" + v"^"'" 
T ^ 2^ 2^ ' 

ay"+i + (l-2a)y" + ay"-i = 

^ - J) (y"+' - 2y"+i + y"-i) + i(y"+i + 2y"+i + y"-i) 



rewrite 


13.41 


in 


the form 










(3.5) 






C^ 


2t 
jD(y"+i 


+ R^^ 

+ 2y"+i + 


- 2y"+i 

T 


+ y"- 
= g", 


-1 



where 



Let 



K = eG - ^C + T (a - j] B. 



^(y" + y"-'), w" = y"-y"-i 
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and rewrite p.5p in the form 

w"+i + w" w"+i - w" 1 

(3.6) C +R + -D(v"+i+y")=g". 

Zt T Z 

Multiplying scalarly both sides of (I3.6|) by 

2(v"+i-v")=w"+i+w", 
we get the equality 

(3.7) ^(C(w"+i +w"),w"+i +w") + i(R(w"+i - w"),w"+i +w") + 

ZT T 

(D(v"+i + v"),v"+i - v") = (g",w"+i +w"). 
Similarly (I2.8p . we have 

(g", w"+i + w") < :^(C(w"+i + w") + ^(^- V", ^")- 
With this in mind, from p.7p it follows 

(3.8) f„+i<£„ + ^(^-V",¥'"), 

where 

f„ = (Dv",v") + i(Rw",w"). 
r 

We formulate the conditions under which the value of Sn determines the square 
of the norm of the difference solution. By virtue of the positivity of operator D it 
is sufficient to require non-negativity of operator R. 

For the energy of operators C and G holds the following coordinate- wise repre- 
sentation 



(Cu, U) = [ A ^ ^ BaUa, ^ BaUa 

Considering 



(Gu,u) 



\ a = l a=l / \a = l / 

p \ ^ 

P Y, [{A-^^^Bo^u^f, IJ^pY {A-^B„u^,B^u^) , 

a— 1 a— 1 

we get 

C <pG. 
Therefore, at cr > 1/4 and 9 > p/2 holds R > 0. We have thus proved the following 
assertion. 

Theorem 3.2. If a > 1/4 and 9 > p/2, than operator R > m H, an additive 
vector scheme 113. 3\) is absolutely stable and for the difference solution holds a priori 
estimate Ii3.8\) with 



^n 



y„^yn-l 



D T 



-hi(R(y"-y"-i),y"-y"-i) 



Proved a priori estimate ( ref (22)) guarantees the stability of the difference 
solution in the half-integer time levels (for v") and is the difference analogue for 
estimate (12.91). 
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4. Generalizations 



We note some of the key research areas that focus on the synthesis and develop- 
ment of the obtained results. 

On the basis of a priori estimate (I3.8|) we obtain the convergence of the solution of 
difference problem p.3p to the solution of differential problem (|1.1|) , (|1.2|) with the 
first order of r. In the standard way [8] we consider the problem for the truncation 
error using a particular scheme for finding the solution at the first time level. 

Instead of p.3p we can use another additive schemes. In the class of vector 
additive schemes, in particular, special attention should be given to the scheme 

/3=1 I3=a+1 

n = 0,l,..., a = l,2,...,p. 

In this case, the time derivative of the several components of the vector solution 
is referred to the upper time-level. Such vector additive schemes are widely used 
[I|[TT| at usual decomposition (|1.4|) . 

Some resources are available when considering more general than ()l.ip . (|1.2|) . 
(|1.5p problems. In our study we restricted ourselves to the simplest problems, where 
operators A, B and the components of splitting of Ba, al, 2, ...,p are constant self- 
adjoint and positive in finite Hilbert space H. These restrictions can be removed 
in some cases, by analogy with the theory of additive schemes for problems (|l.ll) . 
(jl.2p with the usual splitting of (ll.Sp . considering, for example, problems with not 
self-adjoint operators, problem with operator factors [TU1I12) . 

In terms of generalizing the results, the greatest interest is to construct the ad- 
ditive operator-difference schemes for solving the Cauchy problem for evolutionary 
equation (jl.ip in the splitting both operator A and operator B — for the problem 
(jl.ip . (|1.2I) . (|1.4p . (|1.5p . In this case the transition to the new time level is based 
on solving a sequence of problems for equations 

clu 

Ba—j^+AaUa^fa{t), t>0 a = l,2,...,p 

with appropriate initial conditions. 
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